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It  is  well  known  that  a  micro-porous  layer  (MPL)  plays  a  crucial  role  in  the  water  management  of  poly¬ 
mer  electrolyte  fuel  cells  (PEFCs),  and  thereby,  significantly  stabilizes  and  improves  cell  performance.  To 
ascertain  the  exact  roles  of  MPLs,  a  numerical  MPL  model  is  developed  in  this  study  and  incorporated 
with  comprehensive,  multi-dimensional,  multi-phase  fuel-cell  models  that  have  been  devised  earlier. 
The  effects  of  different  porous  properties  and  liquid-entry  pressures  between  an  MPL  and  a  gas  diffusion 
layer  (GDL)  are  examined  via  fully  three-dimensional  numerical  simulations.  First,  when  the  differences 
in  pore  properties  and  wettability  between  the  MPL  and  GDL  are  taken  into  account  but  the  difference  in 
the  entry  pressures  is  ignored,  the  numerical  MPL  model  captures  a  discontinuity  in  liquid  saturation  at 
the  GDL|MPL  interface.  The  simulation  does  not,  however,  capture  the  beneficial  effects  of  an  MPL  on  cell 
performance,  predicting  even  lower  performance  than  in  the  case  of  no  MPL.  On  the  other  hand,  when 
a  high  liquid-entry  pressure  in  an  MPL  is  additionally  considered,  the  numerical  MPL  model  predicts  a 
liquid-free  MPL  and  successfully  demonstrates  the  phenomenon  that  the  high  liquid-entry  pressure  of 
the  MPL  prevents  any  liquid  water  from  entering  the  MPL.  Consequently,  it  is  found  from  the  simulation 
results  that  a  liquid-free  MPL  significantly  enhances  the  back-flow  of  water  across  the  membrane  into  the 
anode,  which,  in  turn,  helps  to  avoid  membrane  dehydration  and  alleviate  the  level  of  GDL  flooding.  As 
a  result,  the  model  successfully  reports  the  beneficial  effects  of  MPLs  on  PEFC  performance  and  predicts 
higher  performance  in  the  presence  of  MPLs  (e.g.,  an  increase  of  67  mV  at  1.5  A  cm-2).  This  study  provides 
a  fundamental  explanation  of  the  function  of  MPLs  and  quantifies  the  influence  of  their  porous  properties 
and  the  liquid-entry  pressure  on  water  transport  and  cell  performance. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

It  has  been  widely  observed  experimentally  that  the  perfor¬ 
mance  of  polymer  electrolyte  fuel  cells  (PEFCs)  is  significantly 
stabilized  and  improved  by  adding  a  micro-porous  layer  (MPL) 
to  the  interface  between  a  gas  diffusion  layer  (GDL)  and  a  cata¬ 
lyst  coated  membrane  (CCM)  [1-8].  Therefore,  the  choice  of  the 
proper  GDL-MPL  combination  is  critical  to  achieve  favorable  PEFC 
performance  and  durability.  An  MPL  is  made  of  a  mixture  of 
carbon  black  powder  and  a  hydrophobic  agent,  usually  polytetraflu- 
oroethylene  (PTFE),  yielding  a  relatively  finer  pore  structure  and  a 
highly  hydrophobic  feature  than  GDLs.  While  the  use  of  MPLs  in 
PEFCs  has  been  a  common  practice  due  to  the  resulting  dramatic 
effects,  the  exact  role  of  an  MPL  is  uncertain  due  to  several  inherent 
difficulties  and,  in  particular,  the  complicated  multi-phase  trans¬ 
port  and  flooding  behaviour  between  the  MPL  and  its  neighboring 
components.  These  components  may  be  either  GDLs  or  catalyst  lay¬ 
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ers  (CLs),  in  of  which  the  morphological  structures  of  the  pores  and 
the  wetting  characteristics  are  all  different  from  those  of  the  MPL. 
Therefore,  recently  published  studies  on  MPLs  focus  mostly  on  a 
more  precise  investigation  and  analysis  of  the  role  of  the  MPL  on 
water  transport  and  the  relevant  PEFC  performance. 

The  effects  of  MPL  composition,  morphology,  and  wettability 
on  reactant  gases  and  product  water  transport  in  PEFCs  have  been 
experimentally  investigated  [9-15].  Wang  et  al.  [9,10]  observed 
the  influence  on  PEFC  performance  of  carbon  powders  that  were 
employed  to  fabricate  MPLs.  They  reported  the  dual  role  of  the 
GDL-MPL  pore  structure  in  the  transport  of  reactant  gases  and 
product  water,  where  the  hydrophobic  pores  in  the  GDL  and  MPL 
facilitate  gas  transport  whereas  the  hydrophilic  pores  allow  a 
pathway  for  liquid  transport.  They  concluded  that  an  appropri¬ 
ate  combination  of  hydrophilic  and  hydrophobic  pores  is  crucial 
for  optimizing  a  GDL-MPL  component  for  the  proper  water  man¬ 
agement  of  PEFCs.  Chen  et  al.  [11]  focused  on  MPL  preparation 
methods,  i.eM  conventional  wet-layer  and  novel  dry-layer  meth¬ 
ods.  Given  the  differences  in  the  surface  morphology  and  pore-size 
distribution  between  GDLs  that  use  either  dry-layer  or  wet-layer 
MPLs,  they  demonstrated  that  the  former  MPLs  exhibit  better  PEFC 
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Nomenclature 

a  water  activity  or  effective  catalyst  area  per  unit  of 

total  volume  (m2  m-3) 

A  area(m2) 

Cp  specific  heat  (kj  kg-1 I<-1 ) 

C  molar  concentration  (mol  m-3 ) 

Dk  mass  diffusivity  of  species,  k  (m2  s-1 ) 

EW  equivalent  weight  of  dry  membrane  (kg  mol-1 ) 

F  Faraday  constant  (96,487  C  mol-1 ) 

i0  exchange  current  density  (A  m-2 ) 

j  transfer  current  density  (A  m-3 ) 

j1  diffusive  mass  flux  of  ith  phase  (kg  m-2  s-1 ) 

J  Leverett  function 

h  enthalpy  per  unit  mass  (kj  kg-1 ) 

k  thermal  conductivity  (W  m  K-1 ) 

kr  relative  permeability 

I<  hydraulic  permeability  (m2 ) 

M  molecular  weight  (kg  mol-1 ) 

m  mass  fraction  (a  scalar  with  no  units) 

n  number  of  electrons  in  electrochemical  reaction  or 

diffusivity  correction  factor 
nc  catalyst  coverage  coefficient 

nd  electro-osmotic  drag  coefficient 

P  pressure  (Pa) 

Pc  capillary  pressure  (Pa) 

RH  relative  humidification  of  inlet 

Ru  universal  gas  constant  (8.314J  mol-1  K-1 ) 
s  stoichiometric  coefficient  in  electrochemical  reac¬ 

tion  or  liquid  saturation 

5  source  term  in  transport  equation 

t  time  (s) 

T  temperature  (K) 

u  fluid  velocity  and  superficial  velocity  in  porous 

medium  (ms-1) 

U0  thermodynamic  equilibrium  potential  (V) 

V  volume  (m3) 

Vceii  cell  potential  (V) 

Greeks  letters 

a  transfer  coefficient 

8t  thickness  of  component  i 

£  volume  fraction  of  gaseous  phase  in  porous  region 

£e  volume  fraction  of  ionomer  phase  in  CL 

y  advection  correction  factor 

A  membrane  water  content  ( mol  H2 0 /mol  SO3  - ) 

A“  relative  mobility  of  phase,  a 

0  phase  potential  (V) 

r\  overpotential  (V) 

6  contact  angle  (°) 

fi  viscosity  (kg  m-1  s-1 ) 

p  density  (kg  m-3) 

pmem  dry  membrane  density  ( kg  m-3 ) 

v  kinematic  viscosity  (m2  s-1 ) 

a  surface  tension  (Nm-1)  or  electronic  conductivity 

(Sm-1) 

r  viscous  shear  stress  (N  m-2) 

k  ionic  conductivity  (S  m-1 ) 

§  stoichiometric  flow  ratio 

Superscripts 

c  cathode  or  capillary 

e  electrolyte 

eff  effective  value  in  porous  region 

mem  membrane 


8  gas 

l  liquid 

ref  reference  value 

s  solid 

sat  saturation  value 

Subscripts 
a  anode 

avg  average  value 

BP  bipolar  plate 

c  cathode  or  capillary 

CL  catalyst  layer 

e  electrolyte 

g  gas  phase 

GC  gas  channel 

GDL  gas  diffusion  layer 

H2  hydrogen 

i  species  index 

in  channel  inlet 

m  mass  equation 

mem  membrane 

N2  nitrogen 

02  oxygen 

ref  reference  value 

t  total 

s  solid 

sat  saturation  value 

w  water 

0  potential  equation 

0  standard  condition,  298.15  K  and  101.3  kPa  (1  atm) 


performance  than  the  latter.  They  hypothesized  that  the  benefi¬ 
cial  effect  of  dry-layer  MPLs  is  due  to  a  relatively  larger  number  of 
meso-pores  that  remain  open  for  gas  transport  for  PEFC  operations 
under  conditions  of  high  humidification.  Ong  et  al.  [12]  examined 
the  effect  of  MPLs,  which  were  based  on  PVDF  (polyvinylidene  fluo¬ 
ride),  as  a  function  of  several  material-dependent  parameters,  such 
as  the  PVDF  concentration,  PVDF/electrically  conductive  filler  ratio, 
and  MPL  thickness.  Ramasamy  et  al.  [13]  emphasized  the  influ¬ 
ence  of  MPLs  on  PEFC  durability.  By  measuring  the  water-retention 
capacities  of  new  and  aged  GDL  samples,  it  was  found  that  GDL 
samples  continuously  lose  their  hydrophobicity  during  long-term 
PEFC  operations  and  that  the  presence  of  an  MPL  mitigates  the 
loss  of  hydrophobicity.  Atiyeh  et  al.  [14]  experimentally  investi¬ 
gated  the  effect  of  MPLs  on  water  transport  and  PEFC  performance. 
The  data  indicated  that  the  presence  of  an  MPL  on  the  cathode 
side  induces  better  overall  PEFC  performance  but  does  not  enhance 
back-diffusion  of  water  from  the  cathode  to  the  anode.  Therefore,  it 
was  concluded  that  the  advantage  of  using  an  MPL  in  PEFCs  is  asso¬ 
ciated  neither  with  an  enhanced  back-diffusion  of  water  nor  with 
a  lowered  net  water  drag  coefficient  through  the  membrane  (this 
coefficient  is  the  moles  of  water  dragged  from  the  anode  through 
the  membrane  to  the  cathode  per  mole  of  protons  transported).  By 
contrast,  with  direct  methanol  fuel  cells  (DMFCs),  Liu  and  Wang 
[15]  demonstrated  that  the  use  of  an  MPL  has  a  dramatic  effect  on 
the  net  water  drag  coefficient  through  the  membrane. 

On  the  other  hand,  the  role  of  MPLs  in  overall  water  transport 
in  PEFCs  was  precisely  investigated  from  a  theoretical  perspec¬ 
tive  [16-19].  Several  MPL  models  have  been  introduced.  Nam  and 
Kaviany  [16]  developed  a  one-dimensional  (ID)  fuel-cell  model  in 
which  an  MPL  was  included  between  the  cathode  CL  and  the  GDL. 
Based  on  capillary  water-transport  theory,  it  was  demonstrated  that 
the  water  transport  through  a  porous  medium  can  be  enhanced  by 
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the  spatial  variation  of  several  pore-structure  parameters  (pore  size, 
hydrophobicity,  etc.).  The  ID  simulation  results  predicted  a  discon¬ 
tinuity  in  the  liquid  saturation  profiles  at  the  interface  between  the 
fine  (MPL)  and  coarse  (GDL)  porous  layers.  It  was  concluded  that 
the  use  of  MPLs  reduces  the  accumulation  of  liquid  water  in  the 
GDL  and  MPL,  thereby  facilitating  water  removal  from  the  cathode 
CL  to  the  GDL.  Pasaogullari  and  Wang  [17]  introduced  a  ID  two- 
phase  PEFC  model  based  on  the  multiphase  mixture  (M2)  model 
previously  developed  by  Wang  and  Cheng  [20].  Based  on  their  ID 
simulation  results,  they  concluded  that  the  presence  of  a  discon¬ 
tinuity  of  liquid  saturation  at  the  interface  between  an  MPL  and 
a  GDL  enhances  the  removal  of  liquid  water  through  the  cathode 
GDL  and  consequently,  reduces  the  liquid  saturation  level  in  the 
cathode  CL.  It  should  be  noted  that  both  the  early  MPL  models 
developed  by  Nam  and  Kaviany  [16]  and  Pasaogullari  and  Wang 
[17]  were  based  on  half-cell  models  and  therefore  were,  unable  to 
capture  water-transport  phenomena  across  the  membrane.  Subse¬ 
quently,  Pasaogullari  and  Wang  [18]  and  Weber  and  Newman  [19] 
developed  ID  full-cell  models  that  accounted  for  both  anode  and 
cathode  GDLs  and  the  membrane.  The  studies  indicated  that  when 
an  MPL  is  employed  between  the  cathode  CL  and  the  GDL,  the  water 
back-flow  across  the  membrane  toward  the  anode  is  enhanced  by 
buildup  of  hydraulic  pressure  in  the  cathode  CL  because  a  much 
higher  capillary  pressure  is  required  to  penetrate  through  the  cath¬ 
ode  MPL.  It  is  reported  that  the  enhanced  water  transport  across 
the  membrane  helps  to  prevent  membrane  dehydration  as  well  as 
serious  flooding  in  the  cathode  GDL. 

Although  significant  efforts  have  been  targeted  at  clarifying 
the  role  of  MPLs  on  overall  water  transport  and  performance,  the 
detailed  function  of  an  MPL  in  increasing  PEFC  performance  and 
durability  is  still  under  debate  and  remains  unclear.  This  can  be 
ascribed  to  a  lack  of  experimental  data  that  clearly  represent  the 
effects  of  MPLs  on  water  transport  in  PEFCs.  Furthermore,  the 
MPL  models  developed  earlier  [16-19]  were  all  based  on  simpli¬ 
fied,  ID,  isothermal  fuel-cell  models  in  which  several  key  effects 
were  missed,  particularly  for  a  non-isothermal  two-phase  situa¬ 
tion.  While  isothermal  two-phase  modeling  can  be  performed  in 
ID,  i.e.,  in  the  through-plane  direction,  as  commonly  done  in  the 
literature  [16-19],  non-isothermal  two-phase  modeling  should  be 
multi-dimensional,  due  to  two  important  additional  phenomena, 
namely,  vapor-phase  diffusion  driven  by  the  thermal  gradient  in  the 
two-phase  region  and  the  heat  pipe  effect  due  to  transfer  of  latent 
heat  from  evaporation  in  the  catalyst  layer  (hot  region)  and  re¬ 
condensation  over  the  land  (cold  region).  Both  phenomena  mainly 
depend  on  the  in-plane  thermal  gradient  between  the  channel  and 
the  land,  which  indicates  that  a  two-dimensional,  cross-sectional 
geometry  considering  both  through-plane  and  in-plane  variations 
is  a  minimum  platform  to  study  non-isothermal,  liquid  water  trans¬ 
port  in  GDL/MPL.  This  study  presents  a  fully  three-dimensional 
(3D),  multi-phase,  fuel-cell  model  that  couples  electrochemical  and 
thermal  phenomena.  In  addition,  numerical  MPL  model  to  assess 


the  effects  of  MPLs  on  water  transport  and  flooding  phenomena  is 
developed  and  incorporated  with  the  fuel-cell  model.  The  objective 
of  this  work  is  threefold: 

( 1 )  to  present  a  numerical  MPL  model  that  can  be  incorporated  with 
existing,  multi-dimensional,  fuel-cell  CFD  (computational  fluid 
dynamics)  models  and  can  enable  a  comprehensive  numerical 
study  of  MPLs; 

(2)  to  investigate  the  effects  of  an  MPL  on  water  and  temperature 
distributions  inside  PEFCs  by  performing  numerical  fuel-cell 
simulations; 

(3)  to  clarify  the  role  of  an  MPL  in  improving  PEFC  performance  and 
provide  a  fundamental  explanation  of  MPL  functions. 

The  paper  is  laid  out  as  follows.  It  starts  with  a  brief  summary  of  a 
multi-dimensional,  two-phase,  fuel-cell  model  that  was  developed 
earlier  [21-23],  along  with  a  detailed  description  and  derivation  of 
the  numerical  MPL  model.  Numerical  simulations  are  carried  out 
and  the  role  of  an  MPL  is  examined  in  detail  under  various  operating 
conditions  of  PEFCs  and  properties  of  GDLs  and  MPLs.  Finally,  major 
conclusions  are  drawn,  that  clearly  elucidate  the  benefits  of  an  MPL 
with  regard  to  water  management  and  PEFC  performance. 

2.  Numerical  model 

2.1.  Multi- dimensional  fuel- cell  model 

The  two-phase,  non-isothermal  PEFC  model  used  in  the  numer¬ 
ical  MPL  study  was  based  on  the  multiphase  mixture  (M2)  model 
developed  by  Wang  and  Cheng  [20].  Tables  1-4  briefly  summa¬ 
rize  the  governing  equations  of  the  fuel-cell  model,  relevant  source 
terms,  electrochemical  correlations,  and  transport  properties  in  the 
electrolyte  phase,  respectively.  Other  details  related  to  the  PEFC 
model  are  available  in  [21-23]. 

2.1.1.  Assumptions  of  model 

Employing  the  M2  formulation  to  describe  two-phase  transport 
in  PEFCs,  the  present  two-phase  PEFC  model  makes  the  following 
assumptions: 

(1)  ideal  gas  mixtures  in  a  single-phase  region; 

(2)  laminar  flow  due  to  small  flow  velocities; 

(3)  isotropic  porous  media  (GDLs,  MPLs,  CLs); 

(4)  two-phase  mist  flow  (i.e.,  homogeneous  flow)  in  GCs,  assuming 
that  channel  flooding  is  minimal  where  tiny  droplets  exist  and 
travel  with  the  gas  velocity  inside  the  GCs. 

2.1.2.  Properties  of  two-phase  transport 

The  two-phase  mixture  properties  are  defined  as  a  function  of 
s  and  (1  -s),  which  denote  the  volume  fractions  of  the  open  pore 


Table  1 

Two-phase  steady-state  PEFC  model:  governing  equations. 


Mass 
Momentum 

Species 

Charge 
Energy 


Governing  equations 

V  ■  (pu)  =  Sm  (1) 

Flow  channels  (Navier-Stokes  equations):  (l/e2)V  ■  ( puu )  =  -Vp  +  V  ■  r  (2) 

Porous  media  (Darcy’s  equations):  pu  =  -(K/vJS/p  (3) 

Flow  channels  and  porous  media:  V  ■  ( Yipmfi )  =  V  ■  [pSDj^Vfm*)]  +  V  ■  [(m^  -  mp/]  +  S,-  (4) 

Water  transport  in  membrane:  V  ■  ((pmem /EW)D™mVX)Mw  -  V  ■  (nd(//F))Mw  +  V  ■  ((I<mem/vl)/VPl)  =  0  (5) 

Proton  transport:  V  ■  (KeffV&e)  +  S<p  =  0  (6) 

Electron  transport:  V  ■  (aeffV<Ps)  -S&  =  0  (7) 

V  ■  ( puCgpT )  =  V  ■  (k&VT)  +  V  ■  ( h°fgplul )  +  ST  (8) 
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Table  2 

Two-phase  steady-state  PEFC  model:  source/sink  terms. 


Source/sink  terms 

Mass 

In  the  CLs:  Sm  =  +  Mw V  ■  (D™em(pmem/EW)VX) 

(9) 

Species 

For  water  in  CLs:  St  =  M,[-V  •  (( nd/F)l )  -  ( sf/nF )] 

(10) 

For  other  species  in  CLs:  Si  =  -Mfsf/nF) 

(11) 

Charge 

In  CLs:  S#  =  j 

(12) 

Energy 

In  CLs:  ST  =j{r)  +  T(dU0/dT ))  + 

(13) 

In  membrane:  Sj  —  l2 /Ke^ 

(14) 

spaces  that  are  occupied  by  the  liquid  and  gas  phases,  respectively 

[20]. 

Density  [20]: 

p  =  pl  -s  +  p8  *(1  -s).  (24) 

In  Eq.  (24),  the  gas  mixture  density,  p8,  as  described  by  the  ideal 
gas  law,  varies  with  the  composition  of  the  mixture  (denoted  here 
by  the  mass  fractions,  {m^}).  That  is 

Velocity  [20]: 
pu  =  plul  +  p£ug. 

Relative  permeability: 
k[  =  s4 


(25) 

(26) 

(27) 


and 

fcf  =  (l-s)4. 

Kinematic  viscosity  [20]: 


(28) 


(29) 


In  Eq.  (29),  v8  is  the  kinematic  viscosity  of  the  gas  mixture  that 
varies  with  the  gas  composition  [24]: 


Xifii 


_ 

~  pg  ~  pgZ^yn  0 

i=i  J  lJ 


where 


^  1  (,  Mi 

“  V8  (  +  MJ 

and 


-1/2 


i+  ^  43 


1/2  (Mp  1/4 

(m, 


/c,  [Nsm  2]  =  < 


f  /ch2  =  0.21  x  10-6T°-66 
fiw  =  0.00584  x  lO^T1  29 
/Cn2  =  0.237  x  10-6T°-76 
.  /Co2  =  0.246  x  10"6r°-78 
Relative  mobility  [20]: 

v‘ 

Xs  =  1  -  X1. 


(30) 


T  in  Kelvin. 


(31) 

(32) 


Table  3 

Electrochemical  correlations. 


EC  Mi  =  chemical  formula  of  species  i 
SiM 2  =  ne~ ,  where  <  s,  =  stoichiometry  coefficient 

l  n  =  number  of  electrons  transferred 


Hydrogen  oxidation  reaction  (HOR)  in  anode  side:  H2  -  2H+  =  2e- 

/  c  \  1/2 

Transfer  current  density  [A  m-3  ] :  j  =  ( 1  -  s)"c ai ^  J  ( 

Surface  overpotential  [V]:  r]  =  (j)s-(j)e 

Oxygen  reduction  reaction  (ORR)  in  cathode  side:  2H20  -  02  -  4H+  =  4e- 

7  c  \ 3/4 

Transper  current  density  [Am-3]:  j  =  -(1  -  sfc ai™{  (  Cq°2  ^  j  exp  fFrj^j 
Surface  overpotential  [V]:  r]  =  0S  -  -  U0 

Thermodynamic  equilibrium  potential:  U0  =  1.23  -  0.9  x  10-3(T -  298.15) 

Temperature  dependence  of  ORR  kinetic  parameter:  airQefc(T )  =  airQefc(353K)  exp  [-jfc  (f  -  35TT5  )  ] 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 
(22) 
(23) 


Table  4 

Transport  properties  in  electrolyte  phase. 
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The  effective  diffusivity  of  species  i  in  the  gas  mixture: 

D?^  =  [S(l -s)]"(T  +  _L)  ,  (33) 

where  Df  and  Df  respectively  denote  the  species  diffusivity  in 
the  gas  mixture  and  the  Knudsen  diffusion  coefficient,  which  are 
expressed  as  follows. 

Species  diffusivity  in  the  gas  mixture  [24]: 

D?  =  1  ~ X|  ,  where 

1  j=n 

j 

_  1.013  x  lO"7  T1-75  /l  1  V/2 

w=  P-(x,1/3  +  xf)2  +  ’ 

Xhw  =  7-07>  Xw  =  12.7,  Xnw  =  17-9-  Xow  =  16.6.  (34) 


The  Knudsen  diffusion  coefficient: 


Df  =  - 


/  8RUT\ 

V  ttM,-  ; 


1/2 


Diffusive  mass  flux  of  the  liquid  phase  [20]: 
jl  =  p/u/  -  Xlpu  =  ^\l\gJVPc, 


(35) 

(36) 


where  the  capillary  pressure,  Pc,  and  the  Leverett  function,  J(s)  are 
computed  as  below: 


pc  =  ps-p'  =  cr  cos  #(£) 


1/2 

J(S), 


(37) 


J(S)  = 


1.417(1  -s)-  2.120(1  -s)2  +  1.263(1  -s)3  if  6C  <  90° 

1.417s- 2.120s2  +  1.263s3  if  9C  >  90° ' 

(38) 


Advection  correction  factor  [20]: 

py)m  |  +  Xgmf) 

1  (sp^J  +  (1  -  s)psrr^) 

2  A3.  The  numerical  MPL  model 

The  development  of  the  numerical  MPL  model  starts  with  the 
balance  of  liquid  pressure  across  the  interface  between  an  MPL  and 
GDL: 

PMPL  =  PGDL •  (40) 

The  overall  liquid  pressure  of  an  MPL  or  GDL  in  Eq.  (40)  is 
separated  into  two  parts,  namely,  the  pressure  of  liquid  water  trans¬ 
portable  by  capillary  action,  Plc,  and  the  liquid-entry  pressure,  Plentry, 
which  represents  a  threshold  pressure  to  penetrate  a  given  porous 
layer  (MPL  or  GDL).  Therefore,  Eq.  (40)  can  be  expressed  in  terms 
of  these  two  pressures  as  follows: 

pi  —  pi  i  pi  —  pi  i  pi  _  pi 

pi  i  pi  —  pi  ’  v  i 

1  c,MPL  ^  1  entry  1  c,GDL 

where  Plentry  represents  the  relative  liquid-entry  pressure  difference 
between  a  GDL  and  MPL  due  to  the  differing  pore  structures  and 
wettabilities. 


Using  the  capillary  pressure  definition  given  by  Eq.  (37),  the  bal¬ 
ance  of  liquid  pressure  can  be  written  in  terms  of  the  capillary  and 
gas  pressures  as  follows: 

-Pc, MPL  +  Fmpl  +  Plentry  =  ~Pc,GDL  +  PqdL’  (42) 

Since  the  gas  pressure  should  be  continuous  at  the  interface 
(P^pL  =  P^DL),  by  using  Eq.  (37),  the  final  form  of  the  liquid  pressure- 
balance  term  in  Eq.  (40)  can  be  expressed  as  two-phase  transport 
parameters  in  terms  of  the  capillary  pressure: 

-Pc, MPL  +  P entry  =  ~^c,GDL 

-<rcos(eMPt)(^)  '  ](slmMPL)  + P'enay  = -a  zos(eCDL)[^)  '  Ksltlt,GDL)  ’ 

(43) 


where  sint<MPL  and  sint  GDL  respectively  denote  the  liquid  saturations 
of  the  MPL  and  GDL  at  the  interface.  Fig.  1,  which  is  a  plot  of  Eq. 
(43)  that  uses  reasonable  MPL  and  GDL  properties,  clearly  illus¬ 
trates  a  discontinuity  in  the  liquid  saturation  across  the  GDL|MPL 
interface.  There  are  two  unknown  variables,  sintiMPL  and  sintGDi,  in 
Eq.  (43);  this  requires  the  solution  of  one  more  equation  that  can  be 
obtained  by  satisfying  the  continuity  of  water  flux  across  the  inter¬ 
face.  The  water-species  equation  for  the  GDL  and  MPL  region  can 
be  simplified  by  neglecting  the  convection  term,  due  to  the  small 
velocity  in  these  porous  layers: 

0  =  V  •  [/ogDf’e-flrV(mf )]  +  V  ■  [(mf  -  m\f],  (44) 

where  the  first  and  second  terms  in  the  right-hand-side  represent 
the  vapor-phase  diffusion  and  capillary  transport  for  the  two-phase 
region.  As  schematically  shown  in  Fig.  2,  the  water  flux  continuity 
at  the  interface  can  be  written  in  a  discretized  form  by  combining 
Eqs.  (37),  (38)  and  (44)  as  follows: 


tnmpi  —  tnw^MPL  —  tnWQDL 


=  PS^,‘ 


Axi/2 


int, MPL  ,  g  /  (&K)  x/xff  </.s  dj 

+  (mgw  -  mlw)- — 2 — A/ A .gcr  cos (6>)  ^ 


), 


Sl  -  Sjnt, MPL 

AXi/2 


AX2/2 


g  l  (£/C)1/2  ,  Sint,GDL-S2 

-  +  (m8  -  mw)  y  XlXga  cos(0)  -  j  ^  ^ 


(45) 


Therefore,  the  values  of  the  liquid  saturation  at  the  interface, 
viz.,  sintiMPL  and  sint<GDL,  can  be  obtained  by  solving  Eqs.  (43)  and 


768 


K.  Kang,  H.Ju  /  Journal  of  Power  Sources  194  (2009)  763-773 


MPL 


GDL 


Fig.  2.  Schematic  of  computational  grid  at  interface  between  the  cathode  MPL  and 
GDL  for  numerical  implementation  of  MPL  model  to  predict  a  discontinuity  in  liquid 
saturation. 


(45).  Finally,  the  additional  water-source  term  for  the  numerical 
MPL  model  has  to  be  applied  to  the  water-species  equation  of  the 
current  two-phase  fuel-cell  model  as  follows: 

f  0  for  single-phase  region 

$mpl  —  T  .  (46) 

l  V  •  mmpi  f°r  two-phase  region 

2.2.  Boundary  conditions  and  numerical  implementation 

The  inlet  velocities  in  the  anode  and  cathode  gas  channels  (GCs) 
can  be  expressed  as  a  function  of  the  anode  and  cathode  stoichio¬ 
metric  ratios  (§a  and  §c,  respectively),  the  PEFC  operating  current 
density  (/),  the  cross-sectional  areas  of  the  anode  and  cathode  GCs 
(Aa  and  Ac ,  respectively),  and  the  densities  of  hydrogen  gas  or  air, 
which  are  respectively  a  function  of  the  anode/cathode  inlet  pres¬ 
sure,  temperature,  and  humidity: 

Ul/2F)AmemMH2  ,  %c(I/4F)AmemMo2 

uin,a  = - — ~A - 1  and  uiri'C  = - — - -•  (47) 

PH2Aa  Po2Ac 

In  addition,  an  isothermal  boundary  condition  (80  °C  for  this 
study)  is  applied  to  the  outer  surfaces  of  the  bipolar  plate  for  tem¬ 
perature  calculations.  The  two-phase  PEFC  model  coupled  with  the 
numerical  MPL  model  described  above  is  numerically  implemented 
with  a  commercially  available,  computational  fluid  dynamics  (CFD) 
package,  namely,  STAR-CD,  through  its  user-subroutine  files  [25]. 


(a) 


Thru -plane 


▼  18pm  30pm 

In-plane  ^  0.5mm  m  210pm  ,12pm  ,12pm  m  180pme  0.5mm 
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Half  land 
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Fig.  3.  Three-dimensional,  single-channel,  PEFC  geometry:  (a)  dimensions  and  (b)  mesh  configuration. 
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Fig.  4.  Comparison  of  simulation  results  in  terms  of:  average  liquid  saturation  (s)  accumulated  in  PEFC  components;  net  water-transport  coefficient  (a)  through  membrane; 
and  cell  performance  (Vcen)  at  1.5  A  cm-2. 


3.  Results  and  discussion 

In  this  study,  a  comprehensive  two-phase,  non-isothermal,  PEFC 
model  is  applied  to  a  typical  fuel-cell  geometry  with  an  MPL  at  the 
cathode  side,  as  schematically  shown  in  Fig.  3(a);  the  mesh  con¬ 
figuration  is  displayed  in  Fig.  3(b).  It  is  seen  in  Fig.  3  that  PEFC 
operation  is  based  on  the  counter-flow  configuration  of  hydrogen 
in  the  anode  GC  and  air  in  the  cathode  GC.  The  properties  of  the 
individual  cell  components  and  the  operating  conditions  are  listed 
in  Table  5  and  the  physical  properties  are  summarized  in  Table  6.  It 
should  be  noted  that  the  two-phase  transport  properties  of  the  cat¬ 
alyst  layer  are  assumed  to  be  same  as  those  of  MPL  to  facilitate  the 
present  MPL  analysis  for  two-phase  transport  phenomena  between 
an  MPL  and  a  GDL. 

To  analyze  the  role  of  an  MPL  in  PEFC  operation,  three  cases 
have  been  defined  for  this  study.  Case  1  eliminates  the  effect  of 
an  MPL  at  the  cathode  side  by  applying  the  properties  of  a  GDL  in 
Table  5  to  the  MPL,  while  the  true  MPL  properties  in  Table  5  are 
employed  for  cases  2  and  3.  The  contrast  between  cases  2  and  3 
is  related  to  the  relative  liquid-entry  pressure  difference  between 


the  GDL  and  the  MPL,  which  is  seen  in  Eq.  (40).  Case  2  consid¬ 
ers  the  difference  in  porous  properties  between  the  GDL  and  the 
MPL  but  neglects  the  effect  of  the  liquid-entry  pressure  difference, 
even  though  a  much  higher  liquid-entry  pressure  is  expected  in  an 
MPL  than  in  a  GDL  due  to  the  finer  pore  structure  of  the  former. 
Therefore,  it  is  expected  in  case  2  that  a  saturation  jump  occurs 
at  the  GDL|MPL  interface  because  the  capillary  pressure  is  con¬ 
tinuous  across  the  interface,  which  means  that  the  liquid  water 
saturation  in  the  MPL  is  much  smaller  than  that  in  the  GDL  because 
of  the  smaller  pore  size  and  greater  hydrophobicity  in  the  MPL. 
On  the  other  hand,  case  3  additionally  considers  a  high  liquid- 
entry  pressure  in  an  MPL,  i.e.,  100  kPa.  Since  the  pore  sizes  of  an 
MPL  are  usually  smaller  than  1  p>m,  it  is  reasonable  to  assume  a 
high  value  of  the  liquid-entry  pressure  of  the  MPL.  Consequently, 
it  is  expected  that  in  case  3,  the  capillary-driven  liquid  water  pres¬ 
sure  in  the  GDL  is  far  below  the  liquid-entry  pressure  of  the  MPL, 
which  prevents  any  liquid  water  from  entering  the  MPL.  These  three 
cases  have  been  simulated  at  two  current  densities,  namely,  0.5  and 
1.5  A  cm-2;  so,  altogether,  six  simulations  are  carried  out  in  this 
study. 
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Fig.  5.  Comparison  of  simulation  results  in  terms  of:  average  liquid  saturation  (s)  accumulated  in  PEFC  components;  net  water-transport  coefficient  (a)  through  membrane; 
and  cell  performance  (Vce//)  at  0.5  A  cm-2. 
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Table  5 

Cell  properties  and  operating  conditions. 


Description 

Value 

Porosity  of  GDL  ( sGdl  ) 

0.7 

Porosity  of  MPL  ( sMpl  ) 

0.6 

Porosity  of  CL(£cl) 

0.6 

Volume  fraction  of  ionomer  in  CLs  (£mc) 

0.23 

Permeability  of  GDL  ( I<Gdl  ) 

1.0  x  10-12  m2 

Permeability  of  MPL  ( I<Mpl  ) 

1.0  x  10-13  m2 

Permeability  of  CL  (I<a) 

1.0  x  10-13  m2 

Hydraulic  permeability  of  membrane  ( Kmem ) 

5.0  x  10-20  m2 

Contact  angle  of  GDL  ( 6Gdl  ) 

110° 

Contact  angle  of  MPL  ( 0Mpl  ) 

120° 

Contact  angle  of  CL  (6Cl) 

120° 

Effective  electronic  conductivity  in  GDL  ( ctGdl ) 

10,000  Sm-1 

Effective  electronic  conductivity  in  MPL  ( crMpi ) 

1000  Sm-1 

Effective  electronic  conductivity  in  CL  (aa) 

1000  Sm-1 

Effective  electronic  conductivity  in  BP  ( crBp ) 

20,0005  m-1 

Thermal  conductivity  of  GDL  (I<gdl) 

5.0  Wm-1  K"1 

Thermal  conductivity  of  MPL  (/<mpl) 

1.0  Wm-1  K-1 

Thermal  conductivity  of  CL  (fcCi) 

1.0  Wm”1  K-1 

Thermal  conductivity  of  BP  ( ksp ) 

20.0  Wm-1  K1 

Thermal  conductivity  of  membrane  (kmem) 

0.95  Wm-1  K-1 

Anode/cathode  inlet  pressure  (Pa/Pc) 

1.5atm/1.5atm 

Anode/cathode  stoichiometry  (£a/£c) 

2/2 

Anode/cathode  inlet  relative  humidity  at  353.15  ( RHa/RHc ) 

100%/100% 

For  cases  1-3  at  1.5  A  cm-2,  Fig.  4  shows  the  variation  of  the 
amount  of  the  average  liquid  water  that  is  accumulated  in  an  indi¬ 
vidual  PEFC  component  along  the  through-plane  direction  with  (a) 
the  cell  performance  and  (b)  the  net  water-transport  coefficient 
through  the  membrane.  The  net  water-transport  coefficient  that  is 
denoted  by  “a”  implies  that  the  net  water  flux  across  the  membrane 
is  normalized  by  the  protonic  flux  as  follows: 

_  nd(//F)  -  ( ( pmem /EW )D™em V A.  +  (Kmem/v')VP') 
a-  l/p  . 

First,  for  case  1  (no  MPL  effect),  the  amount  of  liquid  water 
continuously  decreases  from  the  inside  of  the  cell,  i.e.,  CL,  toward 
the  outside,  i.e.,  GC,  and  no  discontinuity  is  observed  in  the  liquid 
saturation.  On  the  other  hand,  it  is  seen  for  case  2  (MPL  without 
the  liquid-entry  pressure  effect)  that  a  discontinuity  in  the  liq¬ 
uid  saturation  across  the  cathode  MPL|GDL  interface  exists  due 
to  differences  in  the  porous  properties  and  wetting  characteristics 
between  the  MPL  and  the  GDL.  The  average  liquid  saturation  in  the 
MPL  is  lower  than  that  in  the  cathode  GDL  due  to  the  relatively 
higher  hydrophobicity  and  smaller  pore  size  of  the  MPL.  Therefore, 
the  MPL  in  case  2  helps  to  lower  the  amount  of  liquid  water  in 


Table  6 

Physical  properties. 


Description 

Value 

Exchange  current  density  x  Ratio  of  reaction  surface 

1.0  x  109  Am-3 

to  CL  volume  in  anode  side  ( ai JyM 

Exchange  current  density  x  Ratio  of  reaction  surface 

2.0  x  104  Airr3 

to  CL  volume  in  cathode  side  (ai^r) 

Activation  energy  for  the  oxygen  reduction  reaction 

73, 269  J  mol-1 

in  cathode  side  ( Ea ) 

Reference  hydrogen  molar  concentration  (cH2,re/) 

40.88  mol  rrr3 

Reference  oxygen  molar  concentration  (c02)re/) 

40.88  mol  m-3 

Anodic  and  cathodic  transfer  coefficients  for 

aa  =  ac  =  i 

hydrogen  oxidation  reaction  (HOR) 

Cathodic  transfer  coefficient  for  oxygen  reduction 

uc  =  1 

reaction  (ORR) 

Dry  membrane  density  ( pmem ) 

2000  kg  irr3 

Equivalent  weight  of  electrolyte  in  membrane  (EW) 

1.1  kg  mol-1 

Faraday  constant  (F) 

96,487  C  mol”1 

Universal  gas  constant  (Ru) 

8.314  J  mol”1  K-1 

Surface  tension  ( cr ) 

0.0625  Nm-1 

Liquid  water  density  ( pl  (80  °C)) 

972  kg  m-3 

Liquid  water  viscosity  (/ 1 l) 

3.5  x  10-4  Ns irr2 

the  cathode  CL  that,  in  turn,  significantly  reduces  the  back-flow  of 
water  across  the  membrane  from  the  cathode  to  the  anode,  thereby 
rendering  the  net  water-transport  coefficient,  a,  of  case  2  positive 
and  the  highest  among  cases  1-3.  Consequently,  a  larger  portion 
of  excessive  product  water  in  the  cathode  CL  is  transported  toward 
the  cathode  side,  raising  the  amount  of  liquid  water  in  the  cathode 
GDL. 

For  case  3  (MPL  with  the  liquid-entry  pressure  effect),  the  MPL 
is  free  of  liquid  water  (zero  liquid  saturation)  due  to  its  high  liquid- 
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Fig.  6.  (a-c)  Comparison  of  liquid  saturation  contours  in  cross-section  of  anode  GDL 
and  cathode  MPL/GDL  at  fractional  distance  of  0.5  from  cathode  inlet:  /  =  1.5  A  cm-2 . 
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entry  pressure,  where  liquid  water  fails  to  invade  the  MPL  by 
capillary  action  and  thus,  only  water  vapor  exists.  This  indicates 
that  the  MPL  for  case  3  acts  as  a  capillary  barrier,  as  a  result  of 
which  water  transport  through  the  MPL  is  only  driven  by  vapor- 
phase  diffusion  under  the  temperature  gradient.  Consequently,  in 
case  3,  the  buildup  of  liquid  water  in  the  cathode  CL,  viz.,  44%  liquid 
saturation,  is  much  higher  than  in  cases  1  and  2  (25.7%  for  case  1 
and  4.21%  for  case  2),  which  in  turn  significantly  enhances  the  water 
back-flow  across  the  membrane  into  the  anode  for  this  case.  As  a 
result,  for  case  3,  a  is  negative  and  the  lowest  among  cases  1-3.  A 
comparison  of  the  liquid  water  accumulation  and  a  value  in  Fig.  4 
indicates  that  the  amount  of  liquid  water  in  the  anode  side  and  a 
are  the  highest  in  case  3  and  the  lowest  in  case  2. 

It  is  also  seen  in  Fig.  4  that  the  differences  in  liquid  water  profiles 
that  are  predicted  in  cases  1  -3  strongly  influences  cell  performance. 
Case  3  has  the  highest  cell  performance  ofVceu  =  0.54  V  at  1.5  A  cm-2, 
which  is  about  67  mV  higher  than  that  of  case  1  and  115  mV  higher 
than  that  of  case  2.  The  higher  performance  in  case  3  must  be  due 
to  the  greater  back-flow  of  water  across  the  membrane  toward  the 
anode  in  case  3,  which  prevents  dehydration  of  the  anode  side  of  the 
membrane  as  well  as  lowers  the  level  of  GDL  flooding.  Therefore, 
the  comparison  of  cases  1-3  in  Fig.  4  clearly  demonstrates  that  the 
MPL  improves  cell  performance  through  improved  water  manage¬ 
ment  that  is  achieved  by  enhancing  the  back-flow  of  water  across 
the  membrane  into  the  anode.  In  addition,  it  is  also  found  that  the 
saturation  jump  itself  at  the  interface  of  the  MPL  and  the  GDL  does 
not  help  to  improve  PEFC  performance,  which  conflicts  with  the 
conclusions  obtained  from  the  previously  published,  ID,  half-cell 
model  simulations  [16,17]. 

Fig.  5  presents  a  comparison  of  the  average  liquid  water  accumu¬ 
lations  in  the  PEFC  components  at  0.5  A  cm-2.  Although  the  same 
trend  is  observed  as  in  Fig.  4,  the  differences  in  the  cell  performance 
and  liquid  water  accumulation  between  cases  1-3  become  smaller, 
which  agrees  with  the  experimental  observations  where  the  ben¬ 
efits  of  MPLs  are  more  noticeable  with  operation  at  higher  current 
densities.  This  indicates  that  the  use  of  MPLs  helps  to  decrease  the 
ohmic  and  mass-transfer  resistances  inside  a  cell. 


Fig.  6  displays  the  liquid  saturation  contours  in  the  anode  GDL 
and  cathode  MPL/GDL  for  cases  1-3  at  1.5  A  cm-2.  No  discontinu¬ 
ity  is  observed  in  the  liquid  saturation  profile  for  case  1  (Fig.  6(a)) 
because  the  MPL  is  not  employed  in  this  case.  On  the  other  hand, 
when  an  MPL  features  between  the  cathode  CL  and  GDL  but  the  rel¬ 
ative  entry-pressure  effect  is  ignored  (case  2),  the  liquid  saturation 
profile,  as  seen  in  Fig.  6(b),  clearly  shows  that  the  presented  numer¬ 
ical  MPL  model  successfully  captures  a  discontinuity  in  the  liquid 
saturation  contour  at  the  GDL|MPL  interface.  This  discontinuity 
arises  from  the  differences  in  the  porous  properties  and  wettability 
between  the  MPL  and  the  GDL.  As  the  liquid-entry  pressure  effect 
is  further  considered  for  the  MPL  (case  3),  it  is  seen  from  Fig.  6(c) 
that  the  single-phase  MPL  situation,  i.e.,  zero  liquid  saturation  in  the 
MPL,  is  successfully  predicted  by  the  current  numerical  MPL  model. 
Since  the  only  means  for  water  transport  in  the  single-phase  MPL  is 
vapor-phase  diffusion  that  is  driven  by  the  temperature  gradient,  a 
significant  enhancement  of  water  back-flow  across  the  membrane 
into  the  anode  is  expected  in  case  3.  Consequently,  a  comparison 
of  the  liquid  water  saturation  contours  in  cases  1  -3  clearly  shows 
that  the  amount  of  liquid  water  in  the  anode  side  is  highest  in 
case  3. 

Fig.  7  shows  the  water  content  contours  at  1.5  A  cm-2  on  a  CCM 
that  includes  the  anode/cathode  CLs  and  the  membrane.  As  shown 
by  the  expressions  given  in  Table  4,  the  proton  conductivity,  /c, 
increases  with  the  water  content,  k,  in  the  electrolyte  phase  of  the 
CCM  and  thus,  a  lower  content  of  water  is  indicative  of  a  greater 
ohmic  loss  in  PEFCs.  First,  all  three  cases  in  Fig.  7  reveal  that  the 
effect  of  membrane  dehydration  is  more  severe  near  the  anode  side 
of  the  CCM  that  faces  the  channel  region,  which  suggests  a  rela¬ 
tively  lower  value  of  water  content  there.  That  is  because  the  anode 
CCM  near  the  channel  is  prone  to  lose  more  water  toward  the  cath¬ 
ode.  More  importantly,  a  comparison  of  the  water  content  profiles 
of  cases  1-3  clearly  answers  why  the  highest  cell  performance  is 
achieved  in  case  3.  The  presence  of  an  MPL  with  the  liquid-entry 
pressure  effect  (case  3)  considerably  prevents  dry  out  of  the  anode 
CCM  by  enhancing  the  back-flow  of  water  across  the  membrane 
into  the  anode  and  as  a  result,  the  water  content  in  the  anode  side 
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Fig.  7.  (a-c)  Comparison  of  water-content  contours  in  cross-section  of  anode/cathode  CLs  and  the  membrane  at  fractional  distance  of  0.5  from  cathode  inlet:  1=  1.5  Acrrr2. 
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Fig.  8.  (a-c)  Comparison  of  temperature  contours  in  cross-section  of  membrane  electrode  assembly  (MEA)  that  includes  anode/cathode  GDLs,  anode/cathode  CLs,  MPL,  and 
membrane,  at  a  fractional  distance  of  0.5  from  cathode  inlet:  1=  1.5  A  cm-2. 


of  the  CCM  for  case  3  is  much  higher  than  that  for  case  1.  On  the 
other  hand,  it  is  seen  from  a  comparison  between  Fig.  7(a)  and  (b) 
that  the  effect  of  anode  dry  out  is  over-predicted  in  case  2  due  to 
the  under-prediction  of  the  water  flux  across  the  membrane  into  the 
anode.  The  water  content  profiles  in  Fig.  7  verify  the  accuracy  of  the 
present  numerical  MPL  model  and  clearly  demonstrate  the  benefi¬ 
cial  effects  of  an  MPL  on  water  management  and  the  performance 
ofPEFCs. 

Fig.  8  presents  the  temperature  contours  on  the  anode  and  cath¬ 
ode  porous  layers  (GDL  and  MPL)  for  cases  1-3  at  1.5  A  cm-2.  As 
discussed  in  regard  to  Fig.  4,  the  cell  performance  of  case  3  is  higher 
than  those  of  cases  1  and  2  due  to  the  effect  of  the  single-phase  MPL 
in  case  3.  As  a  result,  the  maximum  rise  in  the  temperature  for  case 
3  (3.0  °C)  is  lower  than  those  of  cases  1  and  2  (4.75  °C  for  case  1  and 
5.55  °C  for  case  2).  Therefore,  it  can  be  emphasized  that  the  pres¬ 
ence  of  an  MPL  considerably  affects  the  overall  cell  performance  as 
well  as  the  maximum  temperature  rise  and  liquid  saturation  profile 
inside  a  cell. 

4.  Conclusions 

A  numerical  MPL  model  has  been  developed  and  incorporated 
into  a  comprehensive,  non-isothermal,  two-phase  PEFC  model 
that  was  published  earlier  [21-23].  Three-dimensional,  two-phase, 
non-isothermal  simulations  have  been  performed  to  investigate 
precisely  the  role  of  an  MPL  in  the  distribution  of  temperature  and 
water  as  well  as  cell  performance.  The  following  conclusions  can  be 
drawn. 

(1)  When  the  effect  of  an  MPL  due  to  differences  in  porous  prop¬ 
erties  is  taken  into  account  but  the  effect  of  a  high  liquid-entry 
pressure  is  ignored  (case  2),  the  numerical  MPL  model  accu¬ 
rately  capture  a  discontinuity  in  the  liquid  saturation  at  the 
MPL|GDL  interface.  The  discontinuity  causes  a  decrease  in  the 
amount  of  liquid  water  at  the  cathode  CL  that  in  turn  reduces 
the  water  back-flow  across  the  membrane  from  the  cathode  to 
anode.  Consequently,  the  performance  predicted  in  that  case 
(case  2)  is  even  lower  than  in  the  case  without  the  MPL  (case 


1 )  due  to  the  drier  membrane  as  well  as  the  higher  level  of  GDL 
flooding  in  case  2.  The  result  is  contradictory  compared  with 
experimental  observations. 

(2)  When  a  high  liquid-entry  pressure  (100  kPa)  is  further  applied 
to  the  MPL,  the  numerical  MPL  model  predicts  a  single-phase 
MPL  situation  where  liquid  water  fails  to  enter  the  MPL  due  to 
the  high  liquid-entry  pressure.  Only  water  vapor  is  transported 
via  vapor-phase  diffusion  that  is  driven  by  the  temperature  gra¬ 
dient  through  the  MPL.  The  mechanism  significantly  increases 
the  water  back-flow  across  the  membrane  into  the  anode.  Con¬ 
sequently,  consideration  of  the  liquid-entry  pressure  enables 
the  MPL  model  to  predict  successfully  the  beneficial  effect  of 
MPLs  in  terms  of  PEFC  performance,  where  an  increase  of  67  mV 
in  the  performance  is  achieved  at  1.5  A  cm-2.  The  liquid  satura¬ 
tion  profiles  that  are  predicted  by  the  numerical  simulations 
further  indicate  that  the  MPL  mainly  improves  PEFC  perfor¬ 
mance  by  achieving  a  more  favorable  water  profile  inside  the 
cell.  In  turn,  this  is  realized  by  enhancing  the  water  back-flow 
across  the  membrane  toward  the  anode  side,  avoiding  dehydra¬ 
tion  of  the  membrane,  and  alleviating  the  level  of  GDL  flooding. 

(3)  A  comparison  of  the  numerical  simulation  results  at  0.5  and 
1.5  A  cm-2  clearly  demonstrates  that  the  effect  of  MPLs  is  more 
appreciable  at  the  higher  current  density,  which  accords  with 
the  trend  observed  in  much  experimental  data.  This  indicates 
that  the  beneficial  effect  of  MPLs  on  PEFC  performance  is  due 
to  reduction  in  both  the  ohmic  and  the  mass-transfer  resis¬ 
tance.  The  numerical  MPL  model  successfully  captures  the 
MPL  phenomena  and  demonstrates  the  reduction  in  membrane 
dehydration  and  GDL  flooding  that  follow  from  the  presence  of 
an  MPL. 
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